我国HIV-1流行株CRF07_BC和B’耐药进化研究
【摘要】HIV-1CRF07_BC和HIV-1B亚型是我国艾滋病的主要流行病毒。尽管抗病毒治疗能够延长HIV-1感染者的生存寿命、有效的减少艾滋病的传播,但HIV耐药性的产生不仅导致了抗病毒治疗失败,而且增加了耐药病毒株传播的风险。由于HIV-1在药物和体内免疫的双重压力下产生的耐药性及其病毒进化的高度复杂性,使得现有的病毒耐药进化分析受到一定的局限,特别是针对我国流行的HIV-1非B亚型的研究较为有限。为此,本研究从病毒的整体分子变异进化的角度,首先开发了基于选择压力的相关突变分析软件,对HIV-1CRF07_BC的逆转录酶区相关突变模式、耐药进化通路和我国既往献血HIV-1感染人群病毒基因变异对长期抗病毒治疗的影响进行了系统的分析。该项研究有助于理解HIV-1在抗病毒治疗下特异的进化模式和多态性位点在耐药发生发展中的作用,为改进药物组合的治疗策略提供重要的科学依据。1.基于选择压力的相关突变分析软件的开发相关突变(CorrelatedMutation)是进化生物学中的一个基本概念,在一个可编码蛋白质的基因上,氨基酸突变会受到蛋白质功能的限制。我们开发了一个基于R/Bioconductor的软件包CorMut来发现阳性选择位点之间的相关突变。软件包组合了Ka/Ks比值和相关突变分析。CorMut提供了计算单个位点或特异氨基酸Ka/Ks并发现他们之间相关突变的函数。CorMut提供了三种方法来发现相关突变,包括条件选择压力、互信息和Jaccard系数。计算主要由两个步骤组成：第一,发现阳性选择位点；第二,计算阳性选择位点之间的突变关联。值得注意的是第一个步骤是可选择的。同时,CorMut可以通过构建相关突变网络方便的比较两种条件下的相关突变关系。该软件包已发布在国际的生物信息学平台Bioconductor网站,至今已被他人下载2000余次。(http://www.bioconductor.org/packages/devel/bioc/html/CorMut.html)。2.我国HTV-1CRF07_BC流行株逆转录酶区相关突变模式分析为了解HIV-1CRF07_BC的耐药进化,本研究选择了552例未治疗的HIV-1感染者和261例接受抗病毒治疗的HIV-1感染病人[zidothymidine(AZT)/lamivudine(3TC)/nevirapine(NVP)或AZT/3TC/efavirenz(EFV)],分析了HIV-1逆转录酶区和蛋白酶区突变共变异的改变,采用分层网络的方法展示突变的共变异。首先,识别了逆转录酶区和蛋白酶区的三种类型的特征突变：治疗相关突变、治疗拮抗突变和重叠多态性位点。在6个治疗相关突变(K103N,M184V,Q197K,G190A,Y181C和M230L)和5个重叠多态性位点(A36E,R135I,R277K,L283I和D291E)之间发现了10对显著的相关突变。同时,在逆转录酶区的治疗相关突变(I132L)与与蛋白酶区的重叠多态性位点(L10I)之间发现了一对相关突变。最后,发现逆转录酶区与蛋白酶区的重叠多态性位点在整体上与治疗相关突变之间存在显著的相关性,说明多态性位点对耐药突变位点产生一定的影响。与HIV-1欧美B亚型相比,CRF07_BC表现出独特的突变共变异模式,这种存在于HIV-1CRF07_BC的逆转录酶区和蛋白酶区的一些多态性位点在耐药发生过程中可能发挥着重要作用。该研究对于揭示HIV-1CRF07_BC在抗病毒治疗下特异的进化模式,为今后我国药物组合的改进提供重要的科学依据。3.我国HIV-1CRF07_BC流行株逆转录酶区耐药进化通路分析该部分探讨HIV-1CRF07_BC重组病毒在药物压力下的进化通路,并与B亚型进行比较。选择上述未治疗和治疗的HIV感染者样本病毒的逆转录酶区序列,采用选择压力的方法识别耐药相关位点,使用贝叶斯网络法构建CRF07BC在抗病毒治疗下的耐药进化通路,同时利用国际HIV耐药数据库中的数据构建B亚型在相同治疗方案下耐药进化通路,并与CRF07_BC进行了比较。预测的CRF07_BC主要耐药位点为K103N、Q197K、V179D和Y188L。预测的B亚型主要耐药位点为M184V、K103、Y181C、T69N、G190A,K238T、Y188H和P225H。两者重叠较小,但经典的TAM1(41L,210W和215Y)和TAM2(67N,70R和219E/Q)通路均存在于两个通路,与HIV-1欧美B亚型不同的是预测的CRF07_BC的主要耐药位点中不包含TAM相关突变,并且在核苷类药物相关突变与非核苷类药物相关突变之间存在较强的依赖关系。CRF07_BC形成了以K103N、Q197K、V179D和Y188L为主要耐药位点的独特耐药进化通路,并与B亚型的耐药进化通路相比有较大差异。4.我国既往献血HIV-1感染人群病毒基因变异对长期抗病毒治疗的影响为了研究HIV-1在抗病毒治疗下的耐药进化,尤其是其在长期抗病毒治疗情况下的动态变异以及与抗病毒治疗效果,本研究利用我国最初在河南安徽开展建立的抗病毒治疗回顾性队列,首先选择每个病人接受抗病毒治疗最近时间点和未治疗时的各一条逆转录酶区序列,得到一组病人组成的治疗前后的横断面数据,比较治疗前后的HIV-1逆转录酶区序列,然后识别横断面数据的B’耐药新突变及突变关联性,从而分析了HIV-1的进化特征。从队列的角度探索这些新突变及突变关联的动态变化及其与抗病毒治疗效果的关系。在基于横断面的分析中,我们预测了治疗组的新耐药相关突变位点,共得到42个耐药相关的位点。进一步建立了HIV-1三类治疗特征性突变(治疗相关突变、治疗拮抗突变和重叠多态性位点)在治疗前后的相关突变网络。在治疗相关突变与重叠多态性位点之间发现了17对相关突变,这一结果进一步提示治疗前后阳性选择的重叠多态性位点对于耐药相关位点产生的重要作用。这些相关突变对分为上下两大类,逐渐升高和先升高后降低。研究结果提示HIV-1多态性位点数目与病毒学失败时间呈显著的负相关关系,而多态性位点数目与耐药位点出现的时间负相关并不显著。进一步从队列的角度,观察这些多态性位点是否对与之相关的耐药位点的出现时间产生影响,我们筛选得到了16对多态性位点与治疗相关位点相关突变对对耐药相关位点出现的时间存在显著的影响,进一步说明HIV-1逆转录酶区多态性位点变异对于耐药、病毒学失败的产生一定的影响,这一结果提示在设计抗病毒治疗方案时应考虑HIV多态性位点,从而避免某些耐药性位点的产生,进而选择有效的抗病毒药物。
【作者】李臻鹏；
【导师】马丽英；
【作者基本信息】中国疾病预防控制中心，病原生物学，2014，博士
【关键词】HIV-1；耐药性；耐药模式；进化通路；基因多态性；抗病毒治疗；
【参考文献】
[1]林永.几类约束矩阵方程的迭代解法研究[D].上海大学,基础数学,2014,博士.
[2]张丽群.PCV7免疫抑制实验性哮喘的作用及机制研究[D].重庆医科大学,儿科学,2014,博士.
[3]余朵苟,何世伟.应急物流体系构建研究[J].物流科技,2008,11:1-5.
[4]李广兴.电液伺服负载模拟系统通用软件研究[D].西北工业大学,导航、制导与控制,2004,硕士.
[5]石玉涛.南北地震带南段地震各向异性与介质各向异性特征数值模拟初步研究[D].中国地震局地球物理研究所,固体地球物理学,2014,博士.
[6]李硕征.CNONIX标准符合性测试技术研究与应用[D].北方工业大学,计算机技术,2014,硕士.
[7]温运城.银河巨型计算机工程团队师承研究[D].国防科学技术大学,科学技术哲学,2012,硕士.
[8]卜兴华.基于石墨烯的太赫兹调制效应的研究[D].华中科技大学,光学工程,2013,硕士.
[9]庞鸿飞.伦敦奥运会中美女篮进攻技术对比研究[D].北京体育大学,体育教育训练学,2013,硕士.
[10]韩望.竹材原态多方重组单元与冷压热固工艺研究[D].中国林业科学研究院,木基复合材料科学与工程,2014,博士.
[11]刘其.单晶SiC集群磁流变研磨加工工艺实验研究[D].广东工业大学,机械工程,2014,硕士.
[12]吴楠.欠驱动索杆桁架式机械手爪设计与试验研究[D].哈尔滨工业大学,航空宇航制造工程,2014,硕士.
[13]孔令娜.跨文化视角下民族地区入境旅游者决策行为研究[D].西北师范大学,旅游管理,2013,硕士.
[14]程春林.促进水稻生产社会分工对策研究[D].西南大学,区域经济学,2014,硕士.
[15]李一梅.基于用户体验的手机团购客户端交互界面的优化研究[D].天津大学,工业设计工程,2013,硕士.
[16]李振旺.水下高速航行体机动运动非定常超空泡数值模拟[D].哈尔滨工业大学,航天工程,2013,硕士.
[17]谢桂玲.二语学习者词块运用与语言产出质量的行动研究[D].东北师范大学,外国语言学及应用语言学,2012,硕士.
[18]王麒,毛江鸿,陈佳芸,许斌锋,陆飞,俞凯奇.基于BOTDA的钢结构全尺度应变监测方法及其试验研究[J].传感技术学报,2015,03:443-448.
[19]李凡.潍坊市玄武街工程质量控制研究[D].山东大学,建筑与土木工程,2012,硕士.
[20]达龙.报人林白水平民化新闻思想研究[D].河北大学,新闻学,2014,硕士.
[21]亢剑伟.双PWM变频负载信息前馈控制协调研究[D].河南科技大学,电气工程（专业学位）,2014,硕士.
[22]王伟武.基于演化计算的间断参函数识别方法[D].武汉理工大学,计算机应用技术,2004,硕士.
[23]姚俊,唐广波,刘正东,马宝国.棒线材热连轧过程综合数值模拟系统的开发及应用[J].冶金自动化,2010,05:24-29.
[24]周强强.miR-107通过下调SALL4的表达诱导脑胶质瘤细胞凋亡的研究[D].哈尔滨工业大学,生物学,2013,硕士.
[25]盛鑫.混凝土泵车臂架振动分析与控制[D].吉林大学,固体力学,2014,硕士.
[26]缪迪.农村商业银行农户贷款违约风险研究[D].南京农业大学,金融,2013,硕士.
[27]金宏丽,宋爽,刘佳.北大行施食堂监督员制度实现消费者参与沟通监督[J].高校后勤研究.2007(04)
[28]张朋.写生是一种态度[D].首都师范大学,美术学,2013,硕士.
[29]李贤明.燃烧合成制备FeSi_2、NiTi合金粉体及其表征[D].兰州理工大学,材料学,2014,硕士.
[30]汤威.新农村建设背景下乡镇政府债务问题研究[D].吉林大学,马克思主义发展史,2013,硕士.
[31]石爱丽.河南省秋冬季层状云降水特征及物理机制的研究[D].中国气象科学研究院,大气物理学与大气环境,2004,硕士.
[32]严桥.多种降雨条件下沟坡重力侵蚀过程模型试验[D].大连理工大学,水文学及水资源,2013,硕士.
[33]杨少魁.微光与红外图像彩色融合算法研究[D].中国科学院研究生院（西安光学精密机械研究所）,信号与信息处理,2013,硕士.
[34]李成华.危险货物高温限运政策研究[D].北京交通大学,2015.
[35]何林.海明威作品中的两难模式研究[D].吉林大学,比较文学与世界文学,2004,硕士.
[36]李月.两种膨润土复合颗粒的制备及其处理矿山废水中锰的试验研究[D].辽宁工程技术大学,市政工程,2012,硕士.
[37]张昊.累积复合轧制陶瓷/双相不锈钢复合材料的研究[D].兰州理工大学,材料学,2014,硕士.
[38]于江宁.H集团纵向一体化战略研究[D].西南交通大学,工商管理,2012,硕士.
[39]李新运.环境信息系统建设的若干技术进展[A]..山东环境科学学会2002年度学术论文集[C].:,2003:5.
[40]何茜,李吉跃,沈应柏,陈晓阳,尚富华,胡磊,张志毅.毛白杨杂种无性系叶片δ~(13)C差异与气体交换参数[J].植物生态学报,2010,02:144-150.
[41]杜元伟,孙永河,段万春.主观证据交互式提取及融合方法[J].控制与决策,2011,05:732-736.
[42]郭靖.多核系统下并行实时任务节能调度问题研究[D].东北大学,计算机软件与理论,2010,硕士.
[43]夏元峰.变厚度汽车B柱冲压成形工艺研究及模具设计[D].哈尔滨工业大学,机械制造及其自动化,2013,硕士.
[44]李文凯.开放式基金绩效的评价与实证研究[D].吉林大学,金融,2014,硕士.
[45]彭严严.“第四种权力”与上海参议会制度变迁研究（1946-1949）[D].华东师范大学,政治学理论,2013,硕士.
[46]魏越.基于弹性光网络的网络保护技术的相关研究[D].苏州大学,通信与信息系统,2014,硕士.
[47]木鸿英,赵红,吴巧蓉,胡锐.“三校生”英语教学策略研究[J].昆明冶金高等专科学校学报,2005,02:27-30.
[48]罗永龙.磁控溅射制备Si-Ti负极材料及其性能的研究[D].哈尔滨工业大学,化学工程,2013,硕士.
[49]范锐坚.“中华老字号”咀香园品牌发展战略研究[D].广东财经大学,工商管理（专业学位）,2014,硕士.
[50]顾海东.德国西门子在华企业人力资源的获取与再配置问题研究[D].西南交通大学,工商管理（专业学位）,2013,硕士.
